Incipient Wigner Localization in Circular Quantum Dots 
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We study the development of electron-electron correlations in circular quantum dots as the density is de- 
creased. We consider a wide range of both electron number, N < 20, and electron gas parameter, r s < 18, using 
the diffusion quantum Monte Carlo technique. Features associated with correlation appear to develop very dif- 
ferently in quantum dots than in bulk. The main reason is that translational symmetry is necessarily broken in a 
dot, leading to density modulation and inhomogeneity. Electron-electron interactions act to enhance this modu- 
lation ultimately leading to localization. This process appears to be completely smooth and occurs over a wide 
range of density. Thus there is a broad regime of "incipient" Wigner crystallization in these quantum dots. Our 
specific conclusions are: (i) The density develops sharp rings while the pair density shows both radial and angu- 
lar inhomogeneity. (ii) The spin of the ground state is consistent with Hund's (first) rule throughout our entire 
range of r s for all 4 < N < 20. (iii) The addition energy curve first becomes smoother as interactions strengthen 
- the mesoscopic fluctuations are damped by correlation - and then starts to show features characteristic of the 
classical addition energy, (iv) Localization effects are stronger for a smaller number of electrons, (v) Finally, the 
gap to certain spin excitations becomes small at the strong interaction (large r s ) side of our regime. 
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I. INTRODUCTION 

One of the fundamental quests in condensed matter physics 
research is to understand the effects of strong correlation be- 
tween a system's constituent particles. A simplified system 
of particular interest is the "electron gas", 1 in which conduc- 
tion electrons interact pair-wise via Coulomb forces while the 
effect of atoms is ignored. It is well known that the electron 
gas has a Fermi liquid ground state with extended wave func- 
tions in the limit of high electron density, while when the den- 
sity is decreased, thereby increasing the interaction strength, 
electrons become localized in space and order themselves in 
a "Wigner crystal" phase. 1 The interaction strength is often 
parametrized by the gas parameter r s = {cd/a* B ){l/n) 1 / d 
where n is the electron density, d is the spatial dimension, a* B 
is the effective Bohr radius, and Cd is a dimension-dependent 
constant. 2 In two dimensions, r s — l/agfan) 1 / 2 . The 
physics at intermediate r s continues to offer puzzles, both 
from theory and experiments. There has been numerical ev- 
idence in bulk two 3 -* 4 ^ and thread dimensional (2D and 3D) 
systems that a single transition takes place at r£ ,2D f=s 30-35 
and r^' 3D ~ 100. However, recent work in 2D has predicted 
more complex phases and associated transitions or crossovers 
around these critical values ! 8 ' 9 ' 10 ' 11 ' 12,13 ' 14 While the experi- 
mental evidence is largely inconclusive,-*^ and in particular 
the way in which the transition occurs isn't known experimen- 
tally, the problem has drawn a great deal of attention due to the 
hope of uncovering fundamental aspects of correlation effects. 

Over the past decade or so, small confined systems, such 
as quantum dots (QD), have become very popular for exper- 
imental study ,iZii& Beyond their possible relevance for nan- 
otechnology, they are highly tunable in experiments and in- 
troduce level quantization and quantum interference in a con- 
trolled way. In a finite system, there can not, of course, 



be a true phase transition, but a cross-over between weakly 
and strongly correlated regimes is still expected. There are 
several other fundamental differences between quantum dots 
and bulk systems: (a) Broken translational symmetry in a QD 
reduces the ability of the electrons to delocalize. As a re- 
sult, a Wigner-type cross-over is expected for a smaller value 
of r s . (b) Mesoscopic fluctuations, inherent in any confined 
system, 17,19 lead to a rich interplay with the correlation effects. 
These two added features make strong correlation physics par- 
ticularly interesting in a QD. As clean 2D bulk samples with 
large r s are regularly fabricated these days in semiconductor 
heterostructures, 20 it seems to be just a matter of time before 
these systems are patterned into a QD, thus providing an ex- 
cellent probe of correlation effects. 

Circularly symmetric quantum dots have been a focal point 
of theoretical attention for several years. 21 Early calculations 
using density functional theory (DFT) within the local spin 
density approximation showed a spin density wave (SDW) 
signature for r s as small as 2^^- Unrestricted Hartree-Fock 
(UHF) 2 ! yielded both SDW and charge density wave (CDW) 
features for r s ~ 1. These were initially identified as signa- 
tures of strong correlations related to the analog of a WC in 
a finite system, often called a "Wigner molecule". However, 
SDW and CDW both break the fundamental rotational sym- 
metry of the 2D circular dots. Indeed, later calculations con- 
firmed that these effects are largely artifacts of the approxi- 
mations usedi 23,25 ' 26 ' 27 ' 28 ' 29 ' 30 Projection methods were then 
used to restore symmetries as a second stage of the UHF 
calculations i 31 ' 32 ' 33 ' 34 ,35 DFT has been pushed toward the large 
r s limit through the average spin density approximation. 36 
However, the validity of all these methods - those that reduce 
an interacting problem to an effective noninteracting one - re- 
main questionable, particularly for large r s . 

More demanding computational techniques, which treat the 
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correlations in an exact fashion, have also been applied to this 
problem. For example, exact diagonalization (ED) is accurate 
for small QDs ; 26 ' 30 ' 37 but becomes exponentially intractable 
for dots with more than 6 electrons (N > 6) and r s > 4. Path 
Integral Monte Carlo (PIMC) has also been applied to large r s 
circular dots . 38 i 39 i 40 ' 41 i 42 One study 38 found a crossover from 
Fermi liquid to Wigner molecule behavior at r s w 4, a value 
significantly smaller than the bulk r^ 2D . Another, 40 using dif- 
ferent criteria for the transition, found a two-stage transition 
for r s larger than r^ ,2D . Although PIMC treats interactions ac- 
curately, it has its own systematic and statistical problems; for 
instance, it generates a thermal average of states with differ- 
ent L and S quantum numbers, preserving only S z symmetry. 
Since the energy for low-lying spin excitations becomes very 
small in the low density limit, the constraint on the temper- 
ature in order to access the ground state becomes extremely 
stringent. In the absence of a coherent picture, the role of cor- 
relations in QD remains an open problem. 

The most accurate method for treating the ground state of 
strongly interacting quantum dots, in our opinion, is varia- 
tional Monte Carlo (VMC) followed by diffusion Monte Carlo 
(DMC). This has been carried out in a number of cases for 
high to medium density quantum dots i 43 ' 44 ' 45 ' 46 ' 47 We discuss 
this method in more detail below. Briefly, while there is a sys- 
tematic "fixed-node error", it is considerably smaller than the 
systematic and statistical errors of PIMC. 

In this paper, we present a systematic study of circular 
parabolic QDs over a wide range of interaction strengths us- 
ing an accurate VMC and DMC technique. We demonstrate 
how strong interactions bring out the interesting aspects of 
the Wigner physics in QDs; a short report on some aspects 
appeared in Ref. 48. Our main results are: 

(i) The development of inhomogeneities as the interaction 
strength increases is completely smooth, with no discernible 
special value below r s ss 18 for all N < 20. 

(ii) The density develops sharp rings but remains circularly 
symmetric; the pair densities show both radial and angular in- 
homogeneity. 

(iii) The spin of the ground state follows Hund's (first) rule 
throughout our entire range of r s for all 4 < N < 20: it is 
the maximum consistent with the shell structure, in contrast 
to several previous claims * 23 ' 27,36 - 38 ' 42 We do find many viola- 
tions of Hund's second rule - the value of the orbital angular 
momentum is not necessarily a maximum - but a modified 
second rule holds for small r s . 

(iv) The addition energy curve first becomes smoother as 
interactions strengthen - the mesoscopic fluctuations are 
damped by correlation - and then starts to show features char- 
acteristic of the classical addition energy, indicating incipient 
Wigner crystallization. 

(v) Localization effects are stronger for a smaller number of 
electrons. 

(vi) Finally, the gap to certain spin excitations becomes small 
at the strong interaction side of our regime. 

The organization of the paper is as follows. In Section II, 
we describe the model and parameters for the quantum dots 
we study. Section III discusses our technical tools, VMC 
and DMC. We present our results in Section IV. We focus 



our attention first on the density and pair-densities: it is these 
two quantities that encode rich information on the correlation- 
induced inhomogeneities. We then discuss the energy of the 
ground and excited states. The interesting issue of spin-corre- 
lation in a QD in the large r s limit is addressed at the end of 
the Section. Finally, we present our conclusions in Section V. 



II. MODEL AND PARAMETERS 

We consider quantum dots with N electrons confined 
in a circularly symmetric harmonic potential, T4 on (r) = 
m*oj 2 r 2 /2, using the Hamiltonian 

N / ft 2 \ N 2 i 

where m* is the effective mass of the electrons and e is the 
dielectric constant of the medium. We consider 2D systems, 
so that r 2 = x 2 + y 2 . This is because experimental dots made 
by patterning GaAs/AlGaAs heterostructures have very strong 
confinement in the z-direction so that they are essentially two 
dimensional. The last term in the Hamiltonian is the pairwise 
Coulomb repulsion between electrons. The strength of this in- 
teraction is characterized by the gas parameter, r s (in units of 
the effective Bohr radius a* B ), which is related to the average 
density of electrons, n = J n 2 (r)dr/N, by r s = (7m) -1 / 2 
(in 2D). More physically, r s is essentially the ratio between 
the potential and kinetic energy of the system, justifying its 
identification as the interaction strength. We tune r s by vary- 
ing u> in Vcoii(r); this makes the confining potential more nar- 
row or more shallow, making the average density at fixed N 
larger or smaller, thus controlling r s . 

The confining potential prevents the electrons from flying 
apart from each other, and thus an extra positive background 
charge is unnecessary for the stability of the dot. In the 
limit of weak interaction, the density and number of electrons 
in the dot are related to the strength of the confinement by 
lu = e 2 /(em*rgV / A)- In general, the electron-electron inter- 
actions tend to expand the dot, making the effective co smaller 
than the bare u>. This is quite significant for large r s , and so 
the above simple relation between u> and r s breaks down. 

Circular parabolic confinement is a good description for the 
experimental vertical dots, as well as for the few electron lat- 
eral dots where electrons sit in the central region far from the 
confining gates. 49 In this paper, we assume that the circular 
geometry is preserved even at large r s , leaving the issue of 
irregular dots 47 for future study. 

Throughout this paper we use effective atomic units in 
which the length unit a* B is e/m* times the Bohr radius as, 
and the energy is given in effective Hartrees, H* = m* /e 2 
Hartrees. For GaAs, for example, m* = 0.067m e and 
e = 12 .4, leading to a* B = 98 A and H* = 11.9 me V. 

We have studied the Hamiltonian (Q]) for a wide range of 
parameters. We varied oj between 3.0 and 0.0075 for 2 < 
N < 20. For the range of uj considered, r s lies in the range 
0.4 to 18. 
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III. METHOD 

Circular parabolic dots with noninteracting electrons are 
described by single-particle orbitals, called Fock-Darwin 
(FD) orbitals, specified by principal quantum number n and 
azimuthal quantum number I. The energy of the orbitals is 

E n ,i= {2n+ \l\ + l)tkJ. (2) 

The effective interaction between electrons can be built into 
the single particle problem within the framework of a mean 
field theory, such as the DFT or UHF methods mentioned in 
the introduction. These introduce some interaction effects into 
the single-particle wave functions - for instance, the resulting 
orbitals are more extended than the corresponding noninter- 
acting Fock-Darwin orbitals. 

As a starting point of our calculation, we use Kohn-Sham 
(KS) orbitals (p a (r), from a DFT calculation within the local 
density approximation (LDA). We then construct configura- 
tion state functions (CSFs) that are eigenstates of L, S 2 , and 
S z , each of which is a sum of a product of up- and down-spin 
Slater determinants, 

rn 

<S>f' S (R) Y^ KJ)]!)-. (3) 

3=1 

R = {rx, r 2 , r^v} denotes collective coordinates of the N 
electrons, and D', are Slater determinants of spin-up and 
spin-down KS orbitals, respectively. The coefficients (3ij are 
determined by the requirement that $ i,,s (R) is an eigenstate 
of L, S 2 , andSV 

The many-body trial wave function that we use is 

^' s (K) = jJ2 ^f' S (R) (4) 

1=1 

where the Jastrow factor J — J en J ee J een is the product 
of electron-nucleus, electron-electron and electron-electron- 
nucleus factors. By "nucleus" here we mean the center of the 
harmonic external potential. The detailed form of the Jastrow 
factor can be found in Ref.|50l The independent parameters 
to be optimized are the c, and the parameters in J. The lin- 
ear combination of CSFs builds in the near-degeneracy cor- 
relation in the wave function, whereas the Jastrow factor ef- 
ficiently describes the dynamic correlation that would other- 
wise require a very large number of CSFs. The Jastrow factor 
is so effective that only a small number of CSFs is needed. 

In the very weakly interacting limit, the ground state con- 
sists of simply filling the lowest energy single-particle or- 
bitals. Both the circular and harmonic nature of the external 
potential cause degeneracies in the noninteracting spectrum 
Eq. (f2]l and hence in the noninteracting many-body spectrum 
as well. Thus there is a definite "shell structure" in the ener- 
gies of the many-particle states, much as in atoms. The shells 
are full for N = 2, 6, 12, 20, . . .; for these values of N the 
ground state clearly has L = and 5 = 0. For intermediate 
values of N, the orbitals are filled so that the interaction ef- 
fects yield the lowest energy. For the total spin, Hund's first 



rule, familiar from atomic physics, applies here as well: the 
electrons in the open shell arrange so as to have maximum 
possible spin in order to gain exchange energy. For instance, 
the ground state for N = 9 has S = 3/2. 

For all of the results shown in this paper, we use only CSFs 
constructed from the lowest energy shells consistent with the 
desired L and S: no inter-shell excitations are used. Thus for 
the ground state, the CSFs we include are those for the (pos- 
sibly degenerate) noninteracting ground state. For a closed 
shell and L = S= 0, for instance, there is only one such CSF 
For an open shell, there can be more than one CSF meeting 
our criteria; as an illustration, consider the L = S = state of 
N = 8. Six electrons fill the lowest two shells while two are 
distributed among the third shell's three levels, (n = 1, 1 = 0) 
and (n = 0,1 = ±2). The desired state can be made in two 
ways, by either putting both electrons in the former level or 
putting one in each of the latter (in a singlet state). We always 
include both of these CSFs in our calculations for this state. 
For a given N, we consider all L and S which can be obtained 
with CSFs built from orbitals in the lowest energy shell. 

Selected cases are further checked by including more CSFs. 
For N = 3, 6, 7, 9, and 20, some (L, S) values were studied 
with CSFs which included up to two inter-shell excitations. 
For the range of r s studied here, the energies of these states 
were not significantly changed by including these additional 
CSFs, giving confidence in the accuracy of our results. 

We perform both variational Monte Carlo (VMC) and dif- 
fusion Monte Carlo (DMC) calculations. Both the VMC and 
the DMC energies are upper bounds to the true ground state 
energy, E GS - The VMC energy, £vmc > E GS , is 

/^(R)H^ T (R)dR 
EyMC = /*?,(R)* T (R)dR (5) 

= / dRE L (R)V(R) w J2 ^( R ■ 

^ 8=1 

Here, (a) the Nmc Monte Carlo configurations are sampled 
from P{R) = |*t(R)| 2 / / |*T(R)| 2 dR, and (b) E L (R) = 
*y 1 (R)H* T (R) is the local energy which is constant and 
equal to the true energy in the limit that \&y is an exact eigen- 
state. The variational parameters were optimized using the 
variance-minimization method, 51 and some of the results were 
checked by using two recently developed energy minimization 
methods. 

We use diffusion Monte Carlo (DMC) to project out the best 
estimate of Eqs starting with the optimized ^t- The DMC 
method employs the importance-sampled Green function, 

G(R',R,r) = * r (R')(R'|exp(-WT)|R)/*r(R) , (6) 

to project out 4'o(R)^ , t(R), where \&o(R) is the lowest en- 
ergy state that has the same spatial and spin symmetry as 
^t(R)- As r becomes large, the amplitudes of higher en- 
ergy states decay exponentially compared to that of the ground 
state. However, G(R',R, t) is not known exactly for the 
Hamiltonian Ti of Eq. (Q]i, and a short time approximation 
to exp(— Ht) using the Trotter formula is used repeatedly to 
achieve the desired projection. We use a refined algorithm 
following Ref.|54j which has a very small time step error. 
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The mixed estimator for the DMC energy, 

/ ^(R)HI't(R) dR 
J¥*(R) * r (R)dR 



DMC 



(7) 



equals the ground state energy E . However, mixed estima- 
tors Odmc °f operators that do not commute with the Hamil- 
tonian, such as the density, have errors that are linear in the 
error in ^>t- On the other hand, the extrapolated estima- 
tors 20dmc — Ovmc and Od MC /Ovmc have errors that are 
quadratic in the error in \&y. All of the data shown in this pa- 
per for such operators is based on the extrapolated estimators, 
labeled "QMC", unless explicitly marked "VMC". 

The function which minimizes the energy in the absence of 
any constraints has Bosonic symmetry, but we are interested 
in the lowest Fermionic state. Consequently, compared to the 
Fermionic state of interest, the Bosonic component grows ex- 
ponentially fast. In the absence of statistical noise, one could 
still employ symmetry to obtain the Fermionic energy, but in 
a MC method the estimate for the energy has a statistical error 
that grows exponentially with r. The fixed-node approxima- 
tion prevents this catastrophe by imposing the constraint that 
the nodes of Wq(R) are the same as those of 4 , ^(R). The 
result of this constraint is that the fixed-node DMC energy is 
an upper bound to the true energy. For flexible wave func- 
tions with well-optimized parameters the fixed-node error is 
typically very small. 

In our circular quantum dots, since states with angular mo- 
mentum L are degenerate with those of angular momentum 
—L we are free to construct real wave functions by choosing 
the linear combinations, VP^' 5 + ^^, L ' S and ^!^' S — ^^, L ' S . 
In this case we can use the fixed-node approximation. If in- 
stead we choose to employ states with definite L then the wave 
functions are complex, and we must use the fixed-phase ap- 
proximations^^ the generalization of the fixed-node method 
to complex wave functions. Usually the fixed-phase error is 
comparable to and slightly larger than the fixed-node error. 

Technical considerations limit the present study to approx- 
imately r s < 18. For larger r s , one expects that more CSFs 
should be included in the trial wavefunction because excita- 
tions across the shell gaps produced by the interactions be- 
come more important. For example, for the N — 6 ground 
state (L = S = 0) with lu = 0.01, inclusion of CSFs cor- 
responding to two inter-shell excitations lowers the energy 
from 0.689228(5) to 0.689202(5). However, the variance op- 
timization used for most of our results fails to lower -Evmc 
or £"dmc if there are more than 3-4 CSFs in (though it 
does lower the fluctuations of the energy). For a few cases, 
including moderate and large r s and several N, we have done 
preliminary calculations with higher orbitals by including all 
determinants involving promotion of two electrons across a 
shell gap (e.g. 10 CSFs for N — 20), using two recently devel- 
oped energy optimization procedures3S This, then, allows 
for a change in the nodes of ^ T (R). We find that typically a 
multi-CSF calculation produces a slight decrease in the energy 
for r s > 15, with larger decreases per electron for smaller 
TV. The change in the spatial structure, as in the density and 
pair-density discussed below, is even smaller than that in the 
energy. 



TABLE I: The ground state energy of a circular 2D quantum 
dot obtained by three different computational methods: diffusion 
quantum Monte Carlo (DMC), full configuration interaction (CI), 
and path-integral quantum Monte Carlo (PIMC). N, L, and S 
specify the number of electrons in the dot, their angular momentum, 
and their spin. The energy is given in units of fkv, the characteristic 
energy of the external parabolic confining potential. A = (in 
atomic units) characterizes the strength of the interactions. 
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IV. RESULTS 

Examples of DMC energies from our calculations are pre- 
sented in Table 1. For the purpose of comparison, we also 
include energies obtained using other techniques. The DMC 
energies, that are an upper bound on the true ground state en- 
ergies, are lower than those obtained from the other methods, 
showing the accuracy of the method. 

Further results are organized in six topical areas: the elec- 
tron density, the pair density, real vs. complex wave function, 
the addition energy of the dot, the ordering of the different 
(L, S) states in energy, and the nature of the spin correlations. 
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A. Spatial Density Profile 

The evolution of the radial density of electrons, n(r), as 
r a varies is shown in Fig. Q] The ground state density at 
four different r s is plotted for four values of N (7, 9, 16, and 
20). With increasing r s , the electron-electron repulsion ex- 
pands the system spatially; note that the size of the dot scales 
roughly linearly with r s , as expected given the relation be- 
tween r s and the average density. 

The figure shows that increasing r s causes the density pro- 
file to change dramatically. For small r s , the density is fairly 
smooth, with some weak structure coming from the noninter- 
acting shell structure. We found that the weak structure in 
n(r) for small r s is very similar to that obtained from Fock- 
Darwin orbitals consistent with the shell filling. On the other 
hand, large r s induces strong modulation in n(r), resulting in 
the formation of radial rings in the density profile. The rings 
become sharp with increasing r s . Strikingly, once r s becomes 
larger than 10-12, the number of rings is the same as in the 
classical limit 57 for that N: one ring for N < 5, two rings for 

6 < N < 15, and three rings for larger N's up to 20. Further- 
more, for these larger r s , the average number of electrons in 
the outer ring (obtained by simply integrating the density over 
that region) is mostly consistent with the classical value. 

If one employs wave functions of definite orbital angular 
momentum quantum number L, both the total density and 
the spin densities must be circularly symmetric (in two di- 
mensions). Since our VMC and DMC methods do not break 
this symmetry, we obtain circularly symmetric densities in all 
cases. However, at sufficiently large r s even a small perturba- 
tion of the circular symmetry would be sufficient to pin elec- 
trons; in that case, the density would have the sharp peaks of 
a Wigner localized state. In the absence of a perturbation, the 
signature of localization is manifest in the pair-density, as dis- 
cussed in the next subsection, even when it is absent in the 
density. Alternatively, a symmetry broken density profile is 
obtained when L is not fixed; this case is treated in Section 
IV.C. 

We find that the formation of the rings and the increase 
of their sharpness is completely smooth. To quantify this 
statement, we define a quantity, the "Fractional Peak Height" 
(FPH), that tracks the sharpness of the rings in n(r). Its defi- 
nition and dependence on r s is presented in Fig. [2] The FPH 
is the ratio of the depth of the valley in n(r) between the two 
outer rings to the 'average' height of the rings. The construc- 
tion is shown in the inset of Fig. [2] where the FPH is the 
ratio of the two lengths ABj AC. In Fig. [2] the FPH always 
increases with r s as expected. However, the increase is sur- 
prisingly smooth for the whole range of r s studied. We thus 
infer that its increase is not associated with any special value 
of r s signifying a threshold or 'onset'. The continuous thin 
lines are the best fits of our QMC results to power-law behav- 
ior. Even though the best fit curve is ~ ^/rl for small N (e.g. 

7 and 9), for larger N = 16, 20, it is nearly linear. A large 
FPH signifies stronger radial localization as the rings tend to 
decouple from each other. From Fig. [2] we see that for the 
larger values of r s , small N dots are more strongly radially 
localized for a similar interaction strength. 
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FIG. 1 : Radial electron density n(r) = n-\ (r)+nj (r) in the ground 
state for 4 values of N and four interaction strengths r s . (a) N — 
7, (L, S) = (0,1/2), (b)N = 9, (L,S) = (0,3/2), (c) N = 16, 
(L, S) = (0, 2), and (d) N = 20, (L, S) = (0, 0). As r s increases, 
strong radial modulation develops leading to ring structure in n(r). 
The number of rings for r s > 10 is the same as in the classical limit 
(r s — » oo): 2 for the first two panels and 3 for the last two. [The 
values of u> used are: for N = 7 and 9, lo = 0.8, 0.08, 0.02, and 
0.01; for N = 16, u> = 0.269, 0.06, 0.02, and 0.01; and for N = 20, 
a; = 0.8, 0.04, 0.02, and 0.0075.] 



The FPH is, by construction, some sort of 'peak to valley 
ratio', and the construction used here is not unique. However, 
we have checked for a few cases that a different construction 
(e.g. the ratio between the height of the valley and the outer 
peak) leaves our conclusions unchanged. 
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FIG. 2: Evolution of the sharpness of the rings in n(r) with r s , 
which measures the radial inhomogeneity in the dot. The sharpness 
is quantified in terms of the "Fractional Peak Height" (FPH), the con- 
struction for which is shown in the inset, ¥~PH=AB/AC. The main 
panel shows the r s -dependence of the FPH for the TV = 7, 9, 16, and 
20 ground states (same states as in Fig. 1). The lines show the best 
fit of the data to power-law behavior. The FPH grows very smoothly 
with r s for all TV's without any signature of a special threshold value. 
Note that the rings sharpen faster for smaller N, and that the FPH is 
nearly linear for the larger TV's, 



The smooth and featureless increase of FPH with r s devel- 
ops naturally from the radial structure created by noninteract- 
ing shell effects, which can be physically thought of as Friedel 
oscillations due to radial confinement. Our results show that 
these oscillations smoothly grow into strong inhomogeneities, 
finally leading to radial localization. This suggests that the 
electron-electron interactions are acting on pre-existing oscil- 
lations, namely those caused by noninteracting interference 
effects, and smoothly amplifying them. 

We next turn to explore the role of the electron-electron 
correlations in the density profile. First, we show in Fig. [3] 
the density resulting from four successively better treatments: 
noninteracting, DFT in the LDA approximation, VMC, and 
finally fixed-node QMC (extrapolated estimator). The case 
shown is the N = 20 ground state with r s « 15. As expected, 
we see that the interaction significantly modifies n(r) at this 
large r s , first by greatly expanding the dot at the DFT level 
compared to the noninteracting description, and then by am- 
plifying the inhomogeneous structures (e.g. rings), as seen 
by comparing the DFT and QMC results. We note that the 
difference between the DFT and QMC densities is very large 
compared to the corresponding difference for real atomsr& re- 
flecting the fact that the dot is much more strongly correlated. 
Second, we compare in Fig. [4] the FPH obtained from VMC 
to the full QMC result. Interestingly, the VMC result does 
show a feature at r s w 4, namely a kink in the FPH vs. r s 




FIG. 3: The effect of correlations on the radial density at differ- 
ent stages of approximation for the TV = 20 ground state (L — 0, 
5 = 0) with r a ~ 15 (u> — 0.01). The noninteracting density ufd 
(dashed) from Fock-Darwin orbitals is much too compact and has 
very weak radial structure. In the LDA result, tilda (dotted), the 
mutual electronic repulsion makes the dot greatly expand in the ra- 
dial direction compared to ufd- Both the VMC result (dash-dotted) 
and QMC extrapolated estimate (solid) show much stronger inhomo- 
geneous structure due to correlation that the QMC methods build in 
through the Jastrow Factor. 
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FIG. 4: Growth of the fractional peak height (FPH) for the N = 20 
ground state from the QMC extrapolated estimator (solid symbols) 
compared to the less accurate VMC results (open symbols). Note 
that the VMC estimate alone shows a break point in the slope near 
r a w 4. Thus, a DMC calculation, with its better treatment of the 
interactions, is necessary to produce the smooth behavior of the FPH 
as a function of r s . 



curve, which goes away in the full result. This demonstrates 
the insufficient accuracy of VMC in even the modest r s < 5 
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B. Pair-Density 

The issue of correlation induced localization of the indi- 
vidual electrons can not be addressed by looking at the den- 
sity alone, which is manifestly rotationally symmetric for our 
Hamiltonian. The formation of radial rings in the density im- 
plies, as explained in the previous section, radial localization 
of electrons - a feature special to circularly symmetric con- 
fined systems. This symmetry is, however, broken when one 
of the electrons is held fixed at a particular position. Other 
electrons then organize themselves in the dot so as to mini- 
mize their mutual interaction and kinetic energy; in particular, 
when the repulsion is strong enough, electrons localize at the 
classical positions. Therefore, to address the question of indi- 
vidual electron localization we turn to the pair-density. 

The spin-resolved pair density, g aa i (ro; r), is defined as the 
probability of finding an electron with spin a' at location r 
when an electron with spin a is held fixed at ro, and the total 
pair density is gr = <7ft + .9T1- These quantities detect, in 
addition to radial rings, any angular structure induced by the 
interactions. 
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FIG. 5: Evolution of the spin resolved pair-densities, g acr i (ro; r), 
for the N = 20 ground state (L = S = 0). r s increases from the 
top to the bottom; g^ is on the left while <?ff is on the right; ro is 
chosen near the outer ring of the dot. During the initial increase of 
correlation until r s ~ 5, the 'correlation' and the 'exchange' hole 
become fully developed. A further increase of r B produces short 
range order as seen from the bumps along the outer ring near ro. 
The range of these angular oscillations as well as their amplitude 
increase gradually with r s , suggesting the term "incipient Wigner 
localization" for this regime. 



We present in Fig.[5] gj j and for the N — 20 ground 
state for weak, intermediate, and strong interaction strengths 
(r s 0.4, 5, and 15, respectively). The location ro of the 
up-spin electron is fixed at the outermost local maximum of 
the density, i.e. on the center of the outer ring. The behavior 
of g aa i (ro ; r) for small r s is well known: in this weakly inter- 
acting limit, the 'correlation hole' (the hole around ro in g^) 
is much smaller than the 'exchange hole' (the similar hole in 
gil). This is because an opposite spin electron can come ar- 
bitrary close to the fixed one, while a same spin electron is 
forbidden to be in its vicinity by the Pauli exclusion principle. 
We see that the correlation hole gets bigger with r s , becoming 
comparable in size to the exchange hole by intermediate inter- 
action strength (r s ss 5). At the same time that the correlation 
hole gets bigger, so do the closest peaks in the antiparallel- 
spin pair density but not those in the parallel-spin pair density. 
This is a reflection of the fact that the correlation hole inte- 
grates to while the exchange hole integrates to —1. 

Fig. [5] also shows that increasing r s induces weak angular 
modulation in g aa i (ro; r), in addition to the radial ring struc- 
tures. For small to intermediate r s , this modulation can be 
viewed as Friedel oscillations arising due to the fixed electron; 
this is confirmed by noting that the period of the modulation 
matches the Fermi wavelength of electrons in the outer ring. 
Finally, we find that the angular modulation grows with r s . 
Because of its continuous development toward sharper struc- 
ture, this regime shows "incipient Wigner localization". 

To focus further on this angular modulation, we plot 
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FIG. 6: Evolution of the angular modulation in g aa i (ro; r) along 
the outer ring for the same cases as in Fig. [3] is shown by the 
dotted line and g^ by the dashed line. The solid line represents 
the spin-summed pair-density, gr- The location of the fixed electron 
ro defines 6 — 0. Note the clear Friedel oscillations and the result- 
ing spin density modulation in the middle panel. In the lower panel 
(strongest interactions), a strong feature at 6 = -n is present. 
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9aa'(^o', r) along the outer ring on which the up-spin is fixed 
in Fig. [6] We see that the angular oscillations in gx, g-\\, and 
<7ll are damped and weak in comparison with the radial mod- 
ulation. They grow continuously as a function of r s without 
any threshold value, as in the case of radial modulation. We 
find that the pair-density is almost featureless even for r s sig- 
nificantly greater than 1, while the short-range correlations set 
in for r s > 10. Even at the largest r s studied here, the weak- 
ness of these oscillations suggests that the electrons remain 
essentially delocalized on each ring. 

One of the intriguing features in the spin resolved pair- 
density in Fig. [6] is a bump at 9 = it. At this position, 
which is diametrically opposite to the fixed up-electron, 
decreases while g^ increases compared to their average value 
(the effect of the bump disappears from gr)- We find a sim- 
ilar feature for several TV's with different L, S combinations 
when r s > 5. The feature becomes more pronounced as r s 
increases. We lack a detailed understanding of this "bump" 
at this time; however, we suspect that its origin lies in the na- 
ture of the spin correlation between electrons, which will be 
discussed further in Section IV.F 

We have also studied the pair-density with the position ro 
of the fixed up-electron at different locations. Fig. Q shows 
the spin-resolved pair-densities for the N = 20 ground state 
at r s w 15 with ro chosen on the middle ring (upper panel) 
and at the center 60 of the dot (bottom panel). These plots, 
together with those of Fig. [5] suggest that the angular modu- 
lation is produced primarily in the same ring as the fixed elec- 
tron, while the other rings remain little affected. This is an in- 
dication that even though the FPH at r s s» 15 is substantially 
smaller than its maximum value of 1, the radial localization 
of the electrons is rather strong so that the rings essentially 
decouple from each other. 

We have studied the pair-density for several different N 
over a wide range of r s . As an example, the pair-density for 
a small dot, N — 6, in its ground state (a filled shell case) is 



N=20 




FIG. 7: Spin resolved pair-densities with the fixed electron chosen 
on the middle ring (top panel) or in the center of the dot (bottom 
panel). The parameters are the same as in the bottom panel of Fig. [5] 



shown in Fig. [8]for three values of r s . (The case N — 9, corre- 
sponding to a half-filled shell, was shown in our previous pa- 
per, Ref.|48l) In the classical limit (r a — > oo), two rings are ex- 
pected for N = 6 - a single electron in the center and an outer 
ring containing the remaining five. Note that the r s ss 16.3 
result is quite consistent with this classical structure. The way 
in which the total modulation is shared between and g^ 
near 9 = it is surprising [see Fig. [3d)]: all of the modulation 
is in <7j j while g^ is constant and smaller. We also notice that 
compared to N = 20 for similar r s , the angular oscillations 
in this case are stronger. From our extensive study, we find, 
using results for both pair-density and density (previous Sec- 
tion), that modulation caused by correlation effects is stronger 
when N is small. 



C. Real vs. Complex Trial Wave-function 

In the absence of a magnetic field, one has the choice, as 
discussed in Section|Ill] of working with either a complex or 
real wave-function. (In fact, for a circularly symmetric po- 
tential, one has that choice even in the presence of a mag- 
netic field.) A real is a superposition of the degenerate 
states with angular quantum number L and — L, so is no 
longer an eigenfunction of L. As a result, the angular part of 
the wave-function has sin(m#) and cos(m6') terms instead of 
exp(im9) (m is an integer). While the electron density com- 
ing from a complex wave-function with definite L must be 
circularly symmetric, a real wave-function may have angular 
structure for states with L ^ 0. 

A calculation using real wave-functions which shows angu- 
lar modulation is presented in Fig. [9] The spin-resolved den- 
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FIG. 8: Pair-densities for the N = 6 ground state (L = S = 0, 
a closed shell configuration), gr for a fixed electron on the outer 
ring is shown for (a) r s m 4.5 and (b) r s m 16.3. The corresponding 
angular modulation of the spin-resolved pair-densities is shown in (c) 
and (d). The magnitude of the angular modulation is clearly larger 
than that for iV = 20 at similar values of r s (compare with Fig. [5j: 
electrons are more localized in smaller dots. 
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sity for N — 19 electrons in the L = 3, 5= 1/2 state is shown 
using both real and complex orbitals. We emphasize that the 
strong modulation seen in the top panels is not intrinsically re- 
lated to strong correlation (though we find that the modulation 
amplitude increases with r s ); in particular, the modulation is 
found even for small r s . Furthermore, the modulation disap- 
pears in calculations using complex wave-functions, even for 
large r s , as it must. Oscillations of a similar nature are found 
in the pair density - any modulation in the density must be 
tracked in the pair-density. Thus we conclude that to identify 
the correlation induced inhomogeneities in a L ^ state, it is 
important to use a complex that has a fixed L. 

A disadvantage of using a complex wave-function is, how- 
ever, that the systematic error is typically slightly larger 
(fixed-phase instead of fixed-node approximation). As a re- 
sult, energy estimates are more accurate when *S?t is real. 



D. Addition Energy 

The addition energy, A 2 E(N), is defined as the second dif- 
ference of the ground-state energy with respect to the number 
of electrons, N, on the dot: 

A 2 E(N) = E GS (N + 1) + E GS (N - 1) - 2E GS (N) . (8) 

The addition energy can be accessed experimentally through 
the spacing between conductance peaks in a Coulomb block- 
ade transport measurement ; 17 ! 18 ' 19,21 of the various quantities 
that the we discuss in this paper, A 2 E is the simplest to mea- 
sure experimentally. The leading term in the addition energy 



of quantum dots is the charging energy. Single particle effects 
and corrections to the simple charging model cause the addi- 
tion energy to vary with N. For example, in the noninteracting 
limit of our Hamiltonian (fl}, the spectrum is given by Eq. (0, 
and thus the total many-body energy is a sequence of straight 
line segments of increasing slope. One finds A 2 E(N) = fico 
for the "magic numbers" N = 2, 6, 12, 20, 30, . . . and zero 
otherwise. These special A's correspond to closed shell struc- 
tures for which A 2 E(N) has peaks above its baseline charg- 
ing value because of the extra stability provided by the gap 
between shells. Weak residual interactions beyond charg- 
ing produce further small peaks in A 2 E for half-filled shells, 
N = 4, 9, 16, 25 .... The general occurrence of such peaks 
is referred to as "mesoscopic fluctuations"; they are always 
present in a confined system. 

The addition energy as a function of electron number is 
shown in Fig.[l0] for three values of u>. As N increases at 
fixed u>, the dot grows larger but the confinement potential 
forces the density to increase as well, causing r s to decrease 
slightly with increasing N . In the regime of current experi- 
mental quantum dots (u = 0.28 corresponding to r s w 2), our 
result is similar to previous studies: 21 a base value caused by 
charging modulated by shell effects give strong (weak) peaks 
for filled (half-filled) configurations. Upon increasing the av- 
erage r s (ui — 0.04, middle trace), we find that the strength of 
these peaks weakens, signaling a reduction in the mesoscopic 
fluctuations as a result of electronic correlations. In particu- 
lar, notice that the peak at N = 6 has disappeared completely! 
The fact that the A 2 E curve is smoother for small N indicates 
that correlations are stronger in that regime, a result consistent 
with our inference from densities and pair-densities. 




FIG. 9: Spin-resolved density calculated with a real (top) or complex 
(bottom) wave-function. Up- and down-spin densities for N — 19 
with L = 3, S = 1/2 are on the left and right, respectively. The 
real trial wave-function mixes angular momenta L and — L and so 
breaks the the rotational symmetry, leading to angular modulation in 
the density. The complex trial function has eigenvalue L and shows 
no angular structure (within our statistics). Thus angular modulation 
occurs in the density if the L symmetry is broken by construction; 
such modulation is not due to correlation effects. 
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FIG. 10: Addition energy (normalized by Tkj) as a function of iV 
for three different u and for the classical limit 57 (r s — > oo), As inter- 
actions strengthen because of decreasing u), the mesoscopic fluctua- 
tions in A 2 E become weaker. Note that this happens more readily 
in the small N limit. Features in the ui — 0.01 trace at small iV are 
remarkably similar to those found in the classical limit, showing that 
electrons are nearly localized for small N. 
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Upon further decreasing u, and so increasing r s , A 2 E(N) 
becomes smooth at large N, with a weak remnant of the shell 
effects. However, the behavior for small N has changed com- 
pletely: it is no longer smooth. Note in particular a new peak 
that appears at N = 7. To show the origin of this peak, we 
plot for comparison the addition energy in the classical limit 
using the ground state energy data of Bedanov and Peeters. 57 
Clearly, the nature of fluctuations in the classical limit is very 
different from those in the noninteracting limit: the peak at 
N = 7 can be understood as due to the extra stability of the 
classical configuration. The remarkable similarity of our re- 
sult for small N at r s rs 16 with the classical result strongly 
indicates that electrons in this regime are well localized. 

We emphasize that these changes in the nature of the fluc- 
tuations of the addition energies occur gradually, and happen 
over different r s ranges for small and large N. Our study in- 
dicates that, first, the noninteracting fluctuations die out and, 
then, the classical fluctuations creep in. 



potential splits this degeneracy to favor states with high I. In 
atoms, the gain in exchange energy is insufficient to overcome 
the energy splitting and so orbitals with the same n and I con- 
stitute a shell for the purpose of defining Hund's rules. In 
contrast, for dots, for the range of r s of interest, the exchange 
energy overcomes the shell splitting and so orbitals with the 
same 2n + \l\ define a shell. Since the exchange interaction 
is present in both atoms and dots, agreement with Hund's first 
rule is expected for weakly interacting dots. However, there 
is no reason to expect Hund's second rule to hold, but one 
can enunciate a modified Hund's second rule - states that oc- 
cupy orbitals with high 1 1 1 are favored - which should hold for 
weakly interacting dots. 

Taking N = 9 as an example, we show in Table II results 
for the low-lying states at three values of r s ; the level structure 
and pair densities are shown in Fig.Qj] Note that the ground 
state has 5 = 3/2 as expected from Hund's first rule. At small 
r s , the higher spin states lie at progressively higher energy be- 
cause they involve promotion across one or more shell gaps: 



E. Reordering of States in Energy 

The shell structure present in circular quantum dots in the 
weakly interacting limit suggests, in analogy with atoms, that 
the orbitals within a shell should be filled in the ground state 
so as to follow Hund's rules: the spin should be the high- 
est possible (first rule) and the angular momentum should be 
the maximum consistent with the first rule (second rule). Im- 
plicit in the rules is the definition of a shell, containing orbitals 
that are either exactly or approximately degenerate in energy. 
In atoms, orbitals with the same principal quantum number, 
n, and the same angular quantum number, I, are exactly de- 
generate, while orbitals with the same n but different I have 
an accidental degeneracy in the noninteracting limit because 
of the Coulomb potential. The self-consistent potential splits 
this degeneracy to favor states with low In quantum dots, or- 
bitals with the same radial quantum number, n, and the same 
\l\ are exactly degenerate, while orbitals with the same value 
of 2n+ |Z| have an accidental degeneracy in the noninteracting 
limit because of the harmonic potential. The self-consistent 



TABLE II: The energy (in units of Tkv) of several low-lying states 
(identified by L and S) of a circular 2D quantum dot for N = 9 
and three values of r s . Note the reordering of states in energy as 
the interaction strength increases; in particular, the S — 7/2 state 
becomes nearly degenerate with the ground S= 3/2 state at large r s . 
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FIG. 11: Pair density and level structure for TV = 9 and lo = 0.01 
(r s rs 15.8). The spin-summed pair density, qt, is shown for the 
lowest energy state in each of the five possible spin sectors; a spin- 
up electron is fixed on the outer ring of electrons. The states in the 
top row are the lowest in energy and degenerate within our statistical 
error [(0, 3/2) is the Hund's rule state]. For each of these states, the 
filling of the orbitals in the trial wavefunction is shown. The classical 
configuration for this N is 2 electrons in the center with 7 in the outer 
ring; this configuration is clearly discernible in the pair densities of 
the higher spin states. 
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the kinetic energy cost of such a promotion is too large for any 
interaction effects to overcome. For instance, a transition from 
the (L, 5) = (0, 3/2) state to the (0, 7/2) state involves pro- 
moting the spin-down electrons in the (n, I) = (0, 1), (0, —1) 
orbitals to spin-up electrons in the (n, I) — (0, 3), (0, —3) or- 
bitals. The 5 = 9/2 state requires a further promotion of the 
spin-down electron in (n, I) = (0, 0) to the (1, 1) orbital. 

Note that (L, 5) = (0, 3/2) is the ground state for all three 
values of r s , in agreement with Hund's rules. As expected, 
there are numerous violations of Hund's original second rule 
in the different spin cases but the modified rule discussed 
above holds at small r s . For example, (1) for N = 9 in 
the 5=1/2 sector (Table II), the L = 2 state (with the 
(n, I) — (0, 2) orbital doubly occupied and the (n, I) = (0, —2) 
orbital singly occupied) has a lower energy than the L = 4 
state (with the (n, I) = (0, 2) orbital doubly occupied and the 
(n, I) — (1, 0) orbital singly occupied), , and, (2) for N = 8 
(see Table I), the ground state at r s as 2 is (L, 5) = (0, 1) 
rather than (2, 1). 

As a function of r s , there are two examples in Table II of 
reordering of excitations: First, in the 5=1/2 sector, the L = 
and L = 4 excited states interchange their position by r s as 
6.7, and then by r s as 16 the L = state replaces L = 2 as the 
lowest energy state. Second, as the strength of the interactions 
increases, the 5 = 7/2 excitation becomes progressively of 
lower energy, interchanging with both 5 = 5/2 states by r s as 
6.7 and then with the three 5 = 1/2 states by r s as 16. 

In general, we find that Hund's first rule is very robust: ac- 
cording to our data, it can be used for the ground state spin 
throughout the range 4 < N < 20 and w > 0.01, corresponding 
to r s < 16. For instance, Tables I and II show that the first rule 
works well for N = 6-9, both at small and large r s . In this 
respect we disagree with the PIMC results of Ref.|38 which 
predicted violation for all these cases. We believe the lack of 
sufficient statistical accuracy in Ref. 38 led to that erroneous 
conclusion. Note, for instance, that the statistical error in the 
PIMC energies is often larger than the energy differences be- 
tween different 5 states that we find. 

The only problematic case for Hund's first rule is N = 10. 
Here for uj < 0.04, we find that a state with spin [the state 
(0, 0)] becomes essentially degenerate within the accuracy of 
our calculation with the expected 5 = 1 ground state (2,1). 
Determining the true ground state in the low-density regime 
in this case must await further work. The near degeneracy of 
these two states has been noted previously i 43 i 46 i 61 The reason 
is clear from the level diagram in Fig. QT| the kinetic energy 
difference in moving an electron between the n = 1 and 
orbitals (due to the small splitting between states having dif- 
ferent principal quantum numbers n) very nearly equals the 
exchange energy of one pair of spins. In fact, N = 10 is the 
only N in our range for which promotion of an electron across 
the intra-shell gap results in a gain in the exchange energy of 
just one pair of electrons, explaining why it is the only incon- 
clusive case. 

The single example of a clear violation of the first rule that 
we have encountered is for N — 3 when uo < 0.028: the 
(0, 3/2) state is lower in energy than the Hund's rule ground 
state (1, 1/2). Our results for the energies are essentially the 



same as from a configuration interaction (CI) calculation 37 
and hence not shown here. (For small N, CI calculations pro- 
duce very precise energies.) The first rule violation in this 
case is probably due to the fact that instead of doubly popu- 
lating the spatially localized I = orbital, electrons prefer to 
populate the more delocalized I = 1 orbital. 

We see a trend toward violation of the first rule at our 
largest r s for small N; this may indicate an actual violation 
for r s larger than was possible in this study. For example, for 
N = 9 and u> = 0.01 in Table II, the highly spin-polarized 
state (0, 7/2) becomes degenerate with the usual Hund's rule 
ground state (0, 3/2), to within our numerical accuracy. It is 
clear from the shell structure why 5 = 7/2 is favored: once 
the gain in exchange energy favors promotion of a spin-down 
electron in orbital (0, 1) across two shell gaps to a spin-up 
electron in (0,3), the (0,-1) electron will be promoted to 
(0, —3), thus explaining why 5 = 7/2 becomes lower in en- 
ergy than 5 = 5/2. On the other hand, 5 = 9/2 requires 
promotion across three shell gaps and so is not necessarily fa- 
vored. Similarly, it is expected that for N = 6, the relative 
energy of state (L, 5) = (0, 2) would decrease substantially at 
large r s compared to both the Hund's rule state (0, 0) and the 
fully polarized state (0, 3). For N — 7, the corresponding state 
is (0, 5/2). This is indeed the result seen in Table I. 

Finally, we present the spin-summed pair-densities, §t, for 
different spin states for N — 9 in Fig. QT| It is clear that 
the more polarized states are better localized: as expected, 
exchange acts to keep the electrons apart. This is a rather 
generic feature in our study for a wide parameter range. 



F. Nature of Spin correlation 

The last topic we address is the nature of the spin correla- 
tion in the dots for our larger values of r s . For values of r s by 
which angular modulation in the pair-density has developed, 
the radial localization is already strong, restricting the motion 
of the electrons in the radial direction. Thus, the electronic 
behavior in this limit is perhaps best described in terms of a 
quasi one-dimensional (ID) system on a circular ring. Effec- 
tive spin interactions and the resultant correlations have been 
studied extensively in quasi-lD and 2D over the years i 62 ' 63 ' 64 
Recently, for instance, Matveev and collaborators have shown 
that in a quasi- ID system "ring exchange" processes dominate 
at intermediate r s , leading to novel ground states. 64 

Our results on circularly symmetric quantum dots suggest 
that the the nature of the spin structure depends on the an- 
gular momentum quantum number L. We illustrate this by 
taking the simplest example, a 4-electron dot. Fig. [T2l shows 
the spin-resolved pair-densities for three different low-lying 
states (L, 5) = (0, 1), (2, 0), and (0, 0) for r s - 18. The fixed 
spin-up electron is at 9 = on the single ring. We find that the 
probability of finding the other electrons is maximum at the 
classical locations — 9 = tt/2, tt, and 37r/2. But the spin of 
the electron at these locations depends crucially on the quan- 
tum numbers L and 5 of the state. The ground state is (0, 1) 
(the Hund's rule state), for which the unconstrained electrons 
(2 up and 1 down) are equally likely to occupy any of the re- 
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maining three classical positions. 

More interesting situations occur for the states (2, 0) and 
(0, 0), for which there are an equal number of up and down 
electrons and which differ only in total angular momentum. 
The L = 2 state shows clear antiferromagnetic correlations. 
The L = state is rather unusual: the 9 = tt position is oc- 
cupied by a down electron, while the remaining up and down 
electron are equally distributed over locations 9 = ir/2 and 
3ir/2. We interpret this equal weight of up and down elec- 
trons as representing a spin lying in the plane; thus, this state 
corresponds to a spin-density wave with wave-vector w/2. It 
is interesting to note that the L = state has lower energy 
than the L = 2 state (energies are given in Fig.fT2l) - the anti- 
ferromagnetic state is disfavored. The ordering of these two 
states is opposite to that in the weak r s limit and is another 
example of a violation of Hund's second rule. How the value 
of L is responsible for the spin correlation (for a given S — 0) 
is not settled yet and will be pursued in future studies. Simi- 
lar results were found for other cases of small N as well; for 
large N, there is little angular localization of electrons on each 
ring, making any conclusion unreliable. However, notice that 
the unusual spin correlation here results in an unexpected sur- 
plus of down spins directly opposite the fixed up spin; to that 
extent, the feature here for N — 4 is similar to that at 9 = n 
noted for iV = 20 in Section IV.B. above. 



V. DISCUSSION 

The emerging picture for the correlation-induced inhomo- 
geneities in circular dots appears to be quite different from 




E=0.51682l(6) E=0.51 7267(5) E=0.5 19033(9) 




FIG. 12: Nature of spin-correlation in the three low-lying states for 
N — 4 at large r s (u> = 0.01, r s w 18). Top and middle panels show 
gfl and g^, respectively. Notice that the probability of finding elec- 
trons at a given location depends crucially on the L and S quantum 
numbers. The schematic spin-correlation is given in the lower panel. 
Interestingly, antiferromagnetic correlation occurs for L = 2, S — 
which has the highest energy among these three states. 



that in the bulk. First, the absence of translational symmetry in 
the radial direction introduces radial localization of the elec- 
trons in rings well before individual electrons localize. This 
is clearly against the conventional notion of a single transi- 
tion or cross-over from a weak to strong correlation regime. 
The radial localization can be tracked by the density due to 
the broken translational symmetry; angular localization is re- 
flected in the pair-density - a quantity that does not respect 
the rotational symmetry of the Hamiltonian. It is clear from 
our study that at a given r s , radial localization is stronger than 
angular. This is expected due to the circular geometry. We 
also note here that the rings in the density could in principle 
be directly observed using scanning tunneling microscopy. In 
some of the low-density electron gas systems, resolution may 
be limited by the fact that the electrons are buried below the 
surface. 20 Systems have been developed, however, in which 
the electron gas is near or at the surface , 65 ! 66 ' 67 and it is these 
systems which provide the best opportunities for observation 
of density rings. 

Second, the transition between the weak and strong corre- 
lation regimes is surprisingly broad. In fact, the completely 
smooth evolution of the FPH with r s suggests the absence of 
any cross-over scale at all. Note that the "smoothness" goes 
far beyond the the usual "rounding" of a phase transition in a 
finite system in which a distinct change of slope occurs in the 
cross-over region. 

Third, the correlation strength depends not only on r s but 
also on the number of electrons in the dot: it is typically 
stronger for smaller N, as evident from our results for FPH, 
pair-density, and addition energy. In particular, the nature of 
the mesoscopic fluctuations in addition energy is very differ- 
ent for small and large r s . While they are determined by 
the noninteracting shell-effects for small r s , the structure of 
the classical configurations dictates the fluctuations for large 
r s . As r s increases, first the noninteracting fluctuations are 
smoothed out, and then the classical fluctuations set in. Our 
results indicate that the value of r s where the change takes 
place depends on N. Thus for confined systems, the change 
from the weak to strong correlation regime is strongly non- 
universal. 

Fourth, the ground state spin is consistent with Hund's 
(first) rule throughout the range of our study, 4 < N < 20 and 
r s < 18. At large r s , the excitation energy of certain strongly 
spin-polarized states become small but they do not become 
the ground state. It would be interesting to see if these polar- 
ized states become energetically favorable for even larger r s . 
The extent of inhomogeneity in the dot shows a clear depen- 
dence on the spin state, with stronger localization occurring 
for larger spin polarization. While the degree of electron lo- 
calization is insensitive to the angular momentum L of the 
state (for the small L here), the nature of the spin-correlation 
(its spatial pattern, for instance) is closely connected to L in 
an intriguing manner. 

These conclusions disagree with much of the previous liter- 
ature, the majority of which predicts a single, small cross-over 
scale accompanied by significant deviations from Hund's rule. 
Because our method yields the lowest energies for quantum 
dots to date (the method is strictly variational) and involves 
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less severe approximations, we believe the present results to 
be more accurate. Much of the problem with the earlier results 
stems from either associated approximations or inappropriate 
criteria for the cross-over. 

Finally, we close with a comment on the connection be- 
tween these results for circular dots and those for the bulk 
two-dimensional electron gas. We speculate that with increas- 
ing N, the deep interior of the dot will behave much like the 
bulk system while the radial rings will persist near the bound- 
ary, reflecting the circular confinement. Because of the rings, 
much of the physics discussed here will persist locally near 



the boundary. 
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